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Gebruikte eenheden: gaussische cgs stelsel (zie II.l.4 vor konversie) 


Verwijzingen naar resultaten ín een ander deel beginnen met het nummer van dat 
deel. 


0. Inleiding 


In dít deel worden enige belangrijke versnellingsmechanismen in de radio- 
astrofysica behandeld. Naar hun aard kunnen deze worden onderverdeeld in 
Stochastische Deeltjesversnelling, waarbij een individueel deeltje nu eens 
vertraagd dan weer versneld wordt maar gemiddeld aan energie wint, en 
Reguliere Deeltjesversnelling, waarbij een deeltje voortdurend in energie 
toeneemt. Als voorbeelden van stochastische versnelling behandelen we Fermí- 
versnelling ín hoofdstuk 1, Bêtatronversnelling in hoofdstuk 2 en Versnelling 
door Plasmaturbulentie in hoofdstuk 3. Onder reguliere versnelling wordt 
Schokversnelling in hoofdstuk 4 behandeld en wordt kort ingegaan op 
versnelling door Konstante Elektrische Velden ín hoofdstuk 5. In dit deel 
wordt vooral voortgebouwd op hoofdstuk II.3. 


STOCHASTISCHE DEELTJESVERNELLING 


1. Fermiversnelling 8 


Fermi (1954) heeft voorgesteld dat de kosmische straling in ons melkwegstelsel 
versneld wordt door reflektie aan bewegende interstellaire wolken. Dit 
versnellingsproces is van statistische aard. Per frontale botsing wint een 
deeltje aan energie AE van de orde AE/E » vole, (zie beneden) waarin E de 
deeltjesenergie en Ve de karakteristieke snelheidsdispersie van de bewegende 
“spiegels” is. Omdat er nu VEV VV meer frontale botsingen optreden dan 
“staart”-botsingen (bíj deze laatste verliest het deeltje een energie 
AE/E > v/e), treedt er netto versnelling op (hierin ís v de snelheid van 
het deeltje). De kosmische stralingsdeeltjes zullen proberen in energie- 
equípartitie te geraken met de bewegende wolken die een veel zwaardere massa 
hebben (dus 5 mv? > k zn ‚2 m deelt jesmassa, Mo wolkmassa) . 
Indien dus de deeltjesenergie bij het reflektieproces gemiddeld toeneemt 


volgens 


dE E 
He EE (l.la) 
a 
volgt 
t/t, 
E(t) = E‚ e (1.1b) 


waarin t, de versnellingstijd ís die we onafhankelijk van de deeltjesenergie 
veronderstellen. 

Indien de energetische deeltjes uít het systeem verdwijnen met een 
karakteristieke tíjdschaal ta (die we ook onafhankelijk van de energie 
veronderstellen) wordt de kans dat een deeltje een ouderdom t bereikt gegeven 
door 


=t/t 
P(tjdt =e $ dt/d, (1.2) 


De kans dat een geplt je op moment t een energie E heeft is dan 
(EE) lt, = Aa [E/E) Ì PD 


=t(E)/t 
P(E)dE = P(t(E)) dt(E) = e = dt (E)/t = 
t t /t 
a a dE 
Ee Bt aen 
e ge e 


Het energiespektrum heeft dus de vorm van een machtswet zoals wordt waarge=- 
nomen. Om het waargenomen spektrum van kosmische straling (N(E) a a te 
verklaren moet echter tle, van de orde één zijn; t, wordt bepaald door bijv. 
vernietiging door nucleaire botsingen of ontsnapping uit de melkwegschijf en 
er is geen a priori reden waarom t/t. een dergelijke verhouding zou hebben. 
We zullen zien dat een Fermíiproces werkzaam is bij de versnelling van deeltjes 


in schokken; dit proces staat de laatste jaren sterk in de aandacht. 


Syrovatski heeft een interessante redenering voor het bestaan van een 
universeel deeltjesspektrum bedacht op thermodynamische gronden: Veronderstel 
dat in het versnellingsgebied equipartitie bestaat tussen de energieën in 


snelle deeltjes Wer? turbulente bewegingen Weerd en magnetische velden We, 


zodat dus 
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1 
en = Vb = M= 3 W. B (1.4) 
Onderstel verder dat energie slechts uit het systeem verdwijnt door weglekken 


van snelle deeltjes en dat dit gebeurt op een voldoend grote tijdschaal zó dat 


‘het evenwicht (1.4) steeds gehandhaafd blijft. Indien het totaal aantal deel- 


tjes in de bron No op elk moment gekarakteriseerd kan worden door één 


karakteristieke energie E‚ dan vinden we voor het weglekken met (1.4): 


d(3 N, E) =E d N, 


zodat N, gE en voor het differentiële spektrum van deeltjes dat de bron 
verlaat 
N(E)dE = —dN, = Ee a, 


Het creëren van dit spektrum berust dus op 't voortdurend toevoeren van turbu- 
lente en magnetische energie aan een steeds kleiner wordend totaal aantal 
deeltjes in de bron. Merk op, dat voor een verondersteld evenwicht Wap = 


ÖCH-W op) iepeve (1.4) het uiteindelijke spektrum N(E) a (28) wordt. 


l.l. Gemiddelde energiewinst 


We leiden nu (1.1) eerst af. 

Transformatie van de snelheid vs: y(l,v) is een viervektor zodat in het 
stelsel dat met snelheid v, langs de z-as beweegt de deeltjessnelheid even- 
wijdig aan de z-as gegeven wordt door (zie I1(3.49)) 


î Vv ek NL 
[Lc 
v t mn (1.5) 
Hie 





c2 


en de deeltjesenergie door 
tt = pn 2 í . 
Ct Ue vig de) Ge 


Bij reflektie verandert vj,’ in het meebewegende stelsel van teken Gr 
blijft onveranderd) zodat in het laboratoriumstelsel de snelheid van het 


EA PERS DE En NE PS 


deeltje na de botsing gegeven wordt door (pas (1.5) en (1.6) nogmaals toe 


met Vv, >= Vs Wij >= vi) 





vl + v 2/c?) - 2v 
ws ald c c 
DENS WK ltd 
l - al ia 


1 - 2v v/e? + v‚,2/c? 


Eb El NS A ESE 
ie 1e 4 En zy FET er 2e (1.8) 
ij, ë En 


De verhouding van het aantalkop-(v‚ > ze in (1.7) en (1.8)) tot staart- 
botsingen ís Cv, + v‚) / (v,, - v‚) zodat gemiddeld 


1 + 3 v‚je® 
<y"> a mennen (1.9) 
[1 - zn, Lea - vla - u 2e 

_ N À % _ 4 y Ln 
en Ay = <y"> = y= z (1.10) 
c2 - As ve ec? 


waarbij we bij het tweede gelijkteken Li verwaarloosd hebben t.o.v. vj, . 


Zij \ de gemiddelde afstand tussen de inhomogeniteiten zodat de gemiddel- 
de botsingstijd At = M/v is, dan vinden we uit (1.10) 


bv 2v 
dE elle : (1.11) 
dt cs, 
Vergelijk dít resultaat met (1.1). De “versnellingstijd" ta is dus gegeven 
door 
2 
ben (1.12) 
5 bv? v 
e A 


en is slechts ín het relativistische geval Cv, > Cc) onafhankelijk van de 
energie. 


In het niet-relativistische geval (E = me? + 5 m v?) reduceert (1.11) tot 


(NR) (1.13) 








en ín het extreem-relativistische geval tot 


2 
v 
dE rn: 


ae dak (ER). (1.14) 


De tijd nodig om een deeltje van niet=relativistische snelheid tot relati- 
vistische energie E te versnellen volgt nu uit (1.13) en (1.14) bij benadering 


als 


\ : 
z “Te 7 [Ll + 2n Y] (1.15) 
Cc 


1.2. Spoedhoekverstrooiing 


De elastische botsingen tussen de deeltjes en de bewegende wolken treden op ín 


aanwezigheid van een magnetisch veld. Fermí stelde twee gevallen voor: 





Bij botsingen van het type a treedt reflektie op in gebieden van hogere 
magnetische veldsterkte; het magnetische moment mv, 2/3 ís immers een 
invariante indien de schaal van’ de veldverandering Le >> FG is met 
Ee = v/a, de gyrostraal (Spitzer 1965). Zij a de spoedhoek van het deeltje 
t,o,VvV. het magnetische veld. Door de versnelling neemt MT) voortdurend toe, 
en de spoedhoek dus af. Zodra sin? a < Bral nar zal het deeltje kunnen ont- 
snappen en houdt de versnelling op. Fermiversnelling kan dus slechts optreden 
indien er voortdurend verstrooiing optreedt (met een 
tijdschaal 8. < M/v) zodat de spoedhoek niet systematisch afneemt bij de 
versnelling. Deze verstrooiing zou kunnen worden veroorzaakt door variaties ín 


1 


B met golfgetal k > EG en w > a. zodat p‚2/B niet behouden is. 


Bij botsingen van het type b in gekromde magnetische velden is een 
dergelijk verstrooiend mechanisme niet nodig. 


1,3. Coulombverliezen 


Versnelling kan slechts optreden als de optredende verliezen geringer zijn dan 


if 


de winst. Een belangrijke vorm van energieverlies bij lage energieën zijn de 


„Coulombbotsingen; bij hoge energieën zijn de belangrijkste vormen van energie- 


verlies synchrotronstraling (II-4), inverse Comptonstraling (II.5) en rem- 
straling (II.6). 

Een schatting voor het energieverlies door Coulombbotsingen vinden we als 
volgt: Beschouw een (iehea. relativistisch) deeltje dat een ander (niet- 
relativistisch) deeltje mét massa M en ladíng q ontmoet. Aangezien M aanvanke- 
lijk ongeveer in rust verkeert vinden we de energieoverdracht per botsing uit 

(E‚' At')? 


on OE 


Nu ís het rustraam van het relativistíische deeltje F ' = yq E = ye a/b? en 
At! = b(yv) Ì zodat Ey = e2q2/2Mb2v2 en 


b 
max 2 2 2 a2 

dE  _ nT ef qf n (M) „ — Hef q° n(M) 

dt Ni M vb db Mv AnA (1.16) 
min 


waarin n(M) de deeltjesdichtheid van soort M ís en A de Coulomblogaritme (zie 
II.6.2). De energieverliezen ten gevolge van botsingen met elektronen zijn dus 
veel belangrijker dan die t.g.v. met protonen (protonendichtheid = elektronen” 
dichtheid) zodat tenslotte 


2 2 


dt mv v te 


1.3.1. Intermezzo: Vergelijking van Coulomb-en remstralingsverliezen. 


Het energieverlies ten gevolge van remstraling bíj ontmoeting van een 


(relativistisch) elektron met een proton ís (zie II.,(3./7)) 


6 


22 3 
mes LS (e E‚)? zet ekrarnard (1.18) 


m“ c° b* v 


Het energetisch elektron verliest dus per tijdseenheid 


b 
max 6 6 
HE Kl 3 m2c°b? db = 3 m2 c° b i = 
min 6 en 
spe av 
II(6.5b) 


De remstralingsverliezen (1.19) verhouden zich dus tot de Coulombverliezen 


(1.17) als 
(dE/dt) 2 2 
rem _ ef vj“ 1 
CET or He a (120) 


er zijn dus voor niet-relativistische energieën veel kleiner dan één e?/hc = 
1/137). Voor elektronen met Y X 10% overheersen echter de verliezen door rem- 


straling. 


1.3.2. Injektieenergie 


Uit vergelijking tussen (1.17) en (1,11) volgt dat de Coulombverliezen juist 
bij lage deeltjesenergieën toenemen terwijl de winst t.g.v. Fermiversnelling 
omgekeerd toeneemt voor hoge energieën. Er ís dus een minimale deeltjesenergie 
waarboven Fermi-versnelling werkzaam is: 


’ 


E 5 “lk 
Eos = 400 be A Vo keV (1.21) 
pn ee „2 3 
waarbij we An A = 30 hebben genomen en ke z M(lpc), n= n/(10 ” cm) en 
va = v/(100 km sec !). 
Het Fermimechaníisme kan dus slechts werken als er een voldoende aantal “super- 


thermische" of voorversnelde deeltjes aanwezig is. 


2e Betatronversnellíing 


Een soortgelijk effekt als Fermíiversnelling treedt op in geval van een tijds- 
veranderlijk doch ruimtelijk homogeen magnetisch veld. Indien de magnetische 
veldsterkte periodiek verandert op een tijdschaal die groot is vergeleken bij 
de gyroperiode van de te versnellen deeltjes dan is p‚2/B wederom behouden. De 
hoeveelheid impuls P, zal dus periodiek variëren en een deeltje verandert 
gemiddeld niet van energie. Indien nu echter een verstrooiingsmechanisme 
werkzaam is dat tijdens deze periodieke variaties de deeltjesverdeling 
isotroop houdt zal er wel systematisch magnetische energie in deeltjesenergie 
worden omgezet: 
Beschouw de volgende schematische cyclus 
1. Begin met een isotrope deeltjesverdeling met impuls Po en veldsterkte Bs 


bass Zat, 2 
zodat P‚ 3 Po en p‚, 3 Po” 


2. Versterk het veld tot B, + AB en neem aan dat verstrooiing niet 


belangrijk is; we vinden dan gemiddeld p,° 5 B (1 +2) en 
nn ZE o 
P 73 Po 
3. Houdt het veld op de waarde B, + AB totdat de deeltjesverdeling weer 
isotroop is; dan geldt gemiddeld p 2 = 8 p 2 (3 + 22) en 
bakte ai 7 o 
P// “B Pe ( +3 ). 
d. Laat het veld weer afdemen tot de waarde Bs dan geldt 
2 2AB AB =-1 
pe tg lltg 
5e Houd het veld op de ‘waarde Be tot isotropie bereikt is. Dan geldt dus 


en p‚/ 2 als in 3. 


2 
2 en] EAD enke 
Pi “Po L+g B) Tra earl 
o o 
zodat 5 
AE AB 
E a mt (2.2) 


Coulombinterakties (botsingen) kunnen voor de verstrooiing zorgen indien 
het plasma voldoende dicht en koud ís, maar dan treedt natuurlijk verhitting 
van het gehele plasma op en níet versnelling van een kleine groep. Ook 
magnetische fluktuattes met w > w, kunnen voor de verstrooiing zorg dragen. In 
het bijzonder kan een anisotrope deeltjesverdeling dit soort golven opwekken 
(instabiliteit) en daarbij relaxeren naar een meer isotrope toestand zonder 


noemenswaardig energieverlies. 


3. Plasmaturbulente deelt jesversnelling 


3.l. Plasmaturbulentie 


3.l.l. Inleiding 


Het probleem bij “hydrodynamische turbulentie” is dat altijd de groei- of 


„dempingstijd van individuele wervelingen vergelijkbaar is met de periode 


(rotatietijd). In tegenstelling hiermee is dit vaak niet het geval bij 
“plasmaturbulentie” waar de reële golffrequenties vaak hoog zijn en de groei 
of demping zwak (/y/ << /w/ ; de uitdrukking voor de groeisnelheid van een 
bepaalde golf (k, wy) van soort o wordt gegeven door I1.(3.40). 


In geval van een plasma in thermisch evenwicht met temperatuur T 


=10- 


verwachten we dat een thermisch ruisníveau in de verschillende golfsoorten 


(extra vrijheidsgraden) aanwezig is van de orde 
a 


(27) ° 





KT (3.1) 


waarbij de integraal zich uitstrekt over dat deel van de golf"impuls"ruimte 
waarin Y(«) = 0 is. Voor Langmuirgolven (elektronenplasmagolven) geldt 


ongeveer (zie II(1.4)), K ís de konstante van Boltzmann) 


A 3 
We zKT k (3.2) 
zodat 
ns ï 
a dal ha. | (3.3) 
e D : 


(zie II(1.9)) voor de meeste astrofysische plasma's. (zie betreffende kolom ín 
II.1.3.). 

Vele astrofysische plasma's zijn turbulent, d.w.z. er komen fluctuaties 
in de elektrische — en magnetische velden voor met amplitudes die groter zijn 
dan wat men voor een plasma ín thermisch evenwicht zou verwachten. 

Wanneer de energiedichtheid in deze fluctuaties klein is vergeleken met 
de thermische energie, 

1 wo 


— << 
N, n KT 





<< 1 (3.4) 


dan spreekt men van zwakke turbulentie. In dat geval mag men de turbulente 
veldfluctuaties opvatten als een superpositie van verschillende golven; dus 


voor bijv. de elektrische veldvektor 


E(k,w°(k)) exp (ik. r - tw(k)t) (3.5) 


8E = XL 5 
ok 


waarin elke komponent voldoet aan de lineaire díispersierelatie II(3.22) 
ej 
A (k‚w) = 0. (3.6) 


Wanneer de energiedichtheid van de veldfluctuattes een behoorlijke fraktie 


wordt van de thermische energie van het plasma, w ad 10? n K T, dan spreekt 





dln: rakende danse an “oink dh a” ida 


„il 


men van sterke turbulentie en krijgt men vaak reeds met sterk niet-lineaire 
effekten te maken (bijv. “trapping” van deeltjes in golven); een ontwikkeling 
in golven als in (3.5) ís dan niet geoorloofd en men komt in het regime van de 


"sterke plasmaturbulentie”. 


3.1.2. Random phase benadering 


Een belangrijke, ín principe waarneembare, grootheid in een turbulent plasma 


is de paar=correlatie 


13 (3.7) 


< 8E, (zot) BE,E',E') > 


d 
tussen de componenten van het turbulente electrische veld op plaats r en r' en 
tijdstip t en t'. De haakjes “< >" geven aan dat het hier om een gemiddelde 
grootheid gaat, waarbij men kan denken aan een tijdsgemiddelde over een snelle 
tijdsschaal t = 11E, Altijd geldt < 8E > = 0. 

Wanneer het plasma homogeen is, d.w.z. dat de eigenschappen zoals 
temperatuur en dichtheid niet van de plaats afhangen, en ook geen tijdsaf- 
hankelijkheid vertoont, dan volgt direct dat moet gelden: 


1,3 = LE “Ee, Eet") « (3.8) 


We ontwikkelen de fluctuatie ÒE ín zijn Fourier componenten: 





+o dk i(k.r — wt) 
5E = f ge Á E(k‚w) e (3.9) 
-e (2n)? ’ 


en vinden zo door dit in (3.7) te gebruiken: 


si? EUN ae 
dek dw dk! aw' 
ns < E‚(k,w) E‚(k',w') > x 
2m) (am* ek aa 
ker FWE tk WER) 
x e (3.10) 


Aan deze vergelijking kan alleen voldaan zijn als geldt (Zie II(1.17)) 


<E,(k,w) E(k',w') > = 


mid 


(2) * * ' ' 
Ss An oe E,‚ (k‚w) E, (k‚w) > Sktk!') S(wtw') (3.11) 
TV» TV 
* 
waar we gebruik hebben gemaakt van E‚(-k, u) = Bi (k‚w) (realiteitseis van 


8E). Eigenschap (3.11) van de Fourier componenten van het turbulente 


elektrische veld in een homogeen, tijdsonafhankelijk plasma noemt men de 


random phase approximation. Een speciaal geval van (3.11) is: 


bh 
< E‚(k‚w) E‚(k'‚w') > = lim LEM IE(k,‚w)/?2 Sktk') Elwtw') (3,12) 
T‚Vse 7 V 


* 
De afhankelijkheid van < E‚(k,w) E, (k‚w) > van { k‚w }is niet willekeurig. 
Immers, bij zwakke turbulentie moeten de Fourier componenten toegestane golven 


in het plasma voorstellen, m.a.w. 
E‚(k,w) = EE‚(K) 276 (w- w'(E)) (3.13) 
o 
en (gebruik (3.11) ín (3.10) en II(l.15a)) 


; u) = lin EES (WES (WP 2 Aw - )/V (3.14) 
Vase 5 J 


3.1.3. Fokker-Planck vergelijking 


E, (kw) E 


Indien een hoog nivo aan plasmagolven aanwezig is, dat als zwakke turbulentie 
kan worden opgevat, zullen de resonante deeltjes een voortdurende sterk 
fluktuerende invloed van het golfveld ondergaan. De individuele golfpakket jes 
zijn immers slechts gedurende een korte tijd At (= (aw) 1) en over korte 
afstand Ar (= ca) bj gekorrelleerd (Hierin is Awde frequentiebreedte van 
het golfpakketje en Ak de spreiding ín golfgetallen). 

We zien nu af van andere krachten (dus ook van (konvektieve) bewegingen 
van het plasma) en veronderstellen dat de evolutie van de 
deeltjesverdelingsfunktie ín het turbulente golfveld beschreven kan worden 


door 
E(X t+ Tv) = | E(X - ax; tE) W‚(X — AX; AX) d(AX) (3.15) 


waarin X = (r‚v) de coördinaten in de faseruimte voorstelt en ws AX) de 


overgangswaarschijnlijkheid dat X met AX verandert in tijdsinterval 7. 


ee 


Per definitie geldt 
JW, (Xs AK) d(AK) = 1 


We ontwikkelen (3.15) ín een Taylorreeks en vinden 


òf 





f(Xst) + Tet ……. =S 
fac) tes - a. Ess: ee 
» * AX $ BXAK 1e 
on on 
* [Ws AX) = AK «5 + K(AKAK) «omg veerle (3.16) 
Met de definities 
< AX > z= f d(AX) AXW (Xs AX) , 
ie (3.17) 
< AXAKX > = | d(AX) AXAX Ws A) , 
volgt uit (3.16) 
DE ò AX 92 KAXAD 
rank nd ax £l + Egt CT fl + ee. (3.18) 


de Fokker-Planck vergelijking. 


32e Versnelling door elektrostatische turbulentie 


In het geval van een homogeen plasma zijn de verdelingsfunktie en de coëffi- 
ciënten onafhankelijk van de ruimtelijke coördinaten, zodat (3.18) reduceert 
tot (we gebruiken hier ímpuls i.p.v. snelheid) 








DE à <Ap» 22 <ApAp> 
as haai ir En ks Ke A e Ei „19 
dt (pst) ôp [ T f(p‚t)] + 5 òpòp : [ = f(p,t)] (3 ) 
<Ap» 
Hierin is —— de coëfficiënt van dynamische frictie en <ApApd/ Tt de diffusie- 
coëfficiënt. 


Opgave: Ga na dat de frictiecoëfficiënt inderdaad de verandering beschrijft 
van de gemiddelde impuls van de deeltjesverdeling, en evenzo de diffusie= 


coëfficiënt de verandering van de spreiding van de impulsverdeling. 


=lÂe 


We leiden nu de coëfficiënten af voor elektrostatische turbulenten ín 
afwezigheid van een magneetveld. 

In principe kunnen we <Ap> en <ApAp> berekenen met behulp van de 
bewegingsvergelijking van één deeltje: 

dp 

ar EE (3.20) 
waarbij het electrisch veld 8E (x,t) “gemeten” wordt door het deeltje langs 


zijn ongestoorde baan, d.w.z. 


BE (x‚t) = ôE(x(t), t) 


nn (3.21) 


De impulsverandering wordt als een kleine storing op die baan gezien. Met 
(3.20) volgt direkt: 


t+T tT 
<Ap, Ap =q2 f dt: Jf de” < gr, (zCe'),t') se, (e),e”) > (3.22a) 
t t 


We gaan over op t“ = t' = s als variabele in plaats van t” 


t+T t=t'+T 
<âp, bp > =q f dt’. |J ds <öE, (z(t'),t') GE, (x(t'+e),t'te )D. (3.22b) 
t tet" 


De paarcorrelatie <ör, (z(e'),e') 6E ,(m(t'+te),t'+ts ) heeft een scherp maximum 

voor s = 0, en zal zeer klein worden voor s > To waarbij en de correlatie 
tijd van de turbulente velden is. 

Wanneer de integratie over s zich uitstrekt over een tijdsinterval 

T >> Ts dan kunnen we de integratiegrenzen (t-t!, t=t'+T) zonder meer 


vervangen door (=e, te). Tevens vullen we de Fourier ontwikkeling in: 


t+T ka d3 k dw d° k' dw' 
Ap; Ap = a? | dt' EN mrd mk 
ij t zo (21) * (2x) * 
Xx < E‚(k‚o) E‚(k',u') > 


i {(etk') ox + ik.vt' + ik'.v(t'+s) — iut' — iw'(t'+s) } 
e B: (3.22c) 


=15- 


Door nu gebruik te maken van (3.11) volgt meteen: 


t+T „to d3 k dw ë 

<ApjApd = a? | dt! | ds | E‚(k,W) E‚ (kou) > X 

j t an (2x) * j 
=i(k.v-w)s 

xe (3.22d) 


De integratie over s voeren we uit: 


to =ilk.v-w)s +T =i(k.v-w)s 
ds e = lim Á ds e - 
=e tse -T 
2 sin (wrk.v) T 
= lim MOE 271 ô (w= kev) 


De integratie over t' kan nu ook worden uitgevoerd, en er volgt: 


<Ap, Ap > d° k dw 
ide Ek) EB (kW) > 2 ôlurkew) (3.23) 
(27) * j 


Men ziet dat de Ô funktie die deeltjes selekteert die resonant zijn met de 


w= kev = 0, 


golf: « 


De berekening van <Ap> is iets ingewikkelder. Ten gevolge van de stoor- 
kracht q ÔE (x,t) op een deeltje wijkt dit een beetje af van zijn ongestoorde 
baan x = xo + vt, een stukje ôx. Dan geldt (niet relativistische deeltjes): 


SE, 
EE 


ge 
ct 


waar: 
d 
gr Ra ÔE (zo + ve, t) (3.24) 


2 
en tot in orde lee | volgt dan: 


dp 
gE U ÔE (zot ve + Ôz(t), t) = 
=q { SE (zot vt, t) + öz(t) ‚5 8E (zo + ut, ©} (3.25) 
Met < 8E > = 0 vindt men dan: 
tkT 1 ' e) Kd 1 
<Apd=aqaf de < Öx(t') « 5 E (zo + ve!, EN) D (3.26) 


—16=- 


Een formele oplossing van het stelsel vergelijkingen (3,24) is gegeven door: 
. 
brl!) = AF de" (t'=t") BE (£ + vt", €") (3.27) 


waar we als randvoorwaarden hebben gebruikt: Öx(t) = Ôv(t) = 0 . 


Hiermee wordt (3.26 ): 


2 tt ' ò 
< Apde= a ij dt" ï de"(t'-t") < EEE") «ze GEE!) > (3.28) 


We gebruiken vervolgens: 


5 BE (zot ut', t') =Â BE (x(E") + vlee"), €") = 


Ee: NE «2 
(t'-t") Ôv 


SE (x(t”) + v(t'=t"), £') | (3.29) 
en gebruiken t' — t” = s als variabele i.p.v. t” « 
Dan kunnen we (3.28) schrijven als: 

Capra? Parr Pae ee alt've), toe) Ee 
SE (x(t'), t') D ; (3.30) 
Met dezelfde argumenten als bij de berekening van < Ap, Ap 3 > mogen we de 


integratiegrens (O0, t'=t) voor s uitbreiden tot (O0, ©) . 


De rest van de berekening loopt precies als bij < Ap, Ap, > : 


J 
dì k dw 
ö tkT km % 
< Apde= ik Ni dt' das Sr < EE, w) E(k‚w) > x 
xe 
=i(k.v-u)s 
We schrijven e = cos [(wrk.v)s] — i sin[(wrkev)s]. . Omdat 


* * 
E (-k‚-w) = E (k‚w) is < E (k‚w) E(k‚w) > invariant voor de verwisseling 
(k‚w) > (-k‚-w) , terwijl sin[(w-k.v)s] dan van teken verwisselt. Dít 


betekent dat integratie over d3 k dw de imaginaire term doet verdwijnen. De 


if 


lim sin [(w-k.v)T] 





reële term na integratie over s levert Ees EEE = nT Ô (wek.v) « 
Na integratie over t' volgt dan: TE 
< Ap > 5 d k du 
= 5d gn Sr < Ek) B (k,o) > 18 (uk.v) = 
à < Ap Ap > 
5e nr md (3.32) 


De laatste identiteit volgt uit (3.23). 
Hiermee kunnen we de Fokker-Planck vergelijking (3.19) schrijven als: 





< Ap En > <Ap Ap > 
EN Ed EE pe 
TnT (hop fp) ) +k op op * a £(p) ) (3.33) 
oftewel: 
dE(p) 
òf _ Ò_ 
Be op (°D: Bp (3.34) 
waarin: 
< Ap Ap > 
tk tn sand 
dk dw 
g2 SGD <E (k, WE” (k, w) > n'ê (u-k.v). (3.35) 


Voor elektrostatische golven (met B=o) geldt E° = (E°k) k/k? zodat met 
(3.14) 


dk kk 
a 2 rn g 2 0, es 
D= mia? S Gt < [El 2 > zn 5 WW) - kv) (3.36) 


Vo 


waarin gesommeerd wordt over de golfsoorten co. 
Vergelijking (3-34) met (3,35) beschrijft de stochastische versnelling 


van deeltjes door elektrostatische zwakke turbulentie zonder magneetveld. 
Versnelling door plasmaturbulentie gaat dus gepaard met diffusie in de 


impulsruimte: 


Andain can a indian are as 


rene ie dln era Sie in tt de Sinden 


__=18- 
fp) Wp) 





Pe P Pe jk 
Slechts de resonante deeltjes (w%) = k.v) worden dus versneld en alléên 
voor die deeltjes geldt (3.34)! 


3.2.1. Deeltjesversnelling door Langmuir-turbulentie 


Dispersierelatie voor Langmuirgolven. Volgens II.3.4 moeten we een uitdrukking 


voor 5,3 vinden: 


Beschouw de Vlasovvgl. voor de elektronen in het niet-relativistische geval 


òf e oi 

dt ad Yv. VÉ kad n Ey = 0 j (3.37) 
waarin £Zf (x,v,t) en we voor de versnelling - e E/m hebben genomen 
(elektrostatische golven B (x,t) = 0. We schrijven 
nu f = fo + fj fz + ……eeseseeeese (3.38) 


waarin we de distributiefunktie f ontwikkelen naar machten ín de grootheid E 
die we “klein” veronderstellen. Voor fo kiezen we een Maxwellse verdeling bij 
temperatuur T (dus onafhankelijk van E). We normeren f zó dat 


f dv f(x, v‚ t) = 1 (3.39) 


In eerste orde in E levert (3.37) nu 


Of pa df, 
of na Fouriertransformatie 
E(k‚w) Of, (v)/ Av 
filk, v‚, w) = ie TER (3.41) 


Het verband met de stroomdichtheid wordt gegeven door (ionen vast veronder- 


steld met 1, = Deo =n) 


Az,t) = -en f filx‚v‚t) v d5 v (3,42) 


of in Fourierrepresentatie 


—-19- 


ik,w) = ven ffilk,v,o) vdv. (3.43) 
De geleidingstensor is gedefiniëerd door (zie II(3.16)): 


(k‚w) E‚(k,w) (3.44) 


ijd) = 9e) E, 


Uit (3.41) en (3.43) vinden we dan 


ine? 3 10E Coe, 
jn ee SUL Ger (3.45) 
zodat (zie I1(3.19)) 
v, Òf /òv 
elko) = 8, is Jd EE (3.46) 


Voor de integraal in (3.46) benaderen we als volgt (bedenk dat fo een 


Maxwellverdeling is en we zijn op zoek naar golven met w/k >> Va) 


v‚ of /òv, Bi. v L 
je EI = k.v/w 5 „fare ® =k.v/w) je 
k.v k.v 2 v‚k kv 
sels jen eG ge (+2) = 
k2 vv 2 k k, v 
=- [ô eeN 
ij oe w2 


Voor het gemak kiezen we k// ? en vinden dan voor II.3.18: 


5 ne kv k2c? 
ii doge E- tT (3.47a) 
w_2 Ik dv, 2 
A33 = 1 -— er [1 kare me (3.47b) 
en overigens nul. Hierin 
2 = 2 A3 
ten, Seft red (3.48) 


De mogelijke oplossingen worden nu gegeven door det A,, 2 A= 0; er zijn dus 


ij 


twee oplossingen met dispersierelatie 


) = oe? + k2e2)*, (3.49) 


=20- 


de elektromagnetische golven, en één oplosisng 
UE) mla rale DS (3.50) 
en pe te ú 
de zogenoemde Langmuirgolven (elektrostatische golven). 
Rl). We leiden nu de verhouding tussen totale en elektrische 


golfenergiedichtheíid voor Langmuirgolven af, gegeven door Rr) (zie 
II(3.41)). 











We vinden 

„A is 

WB 2 = 2 A11 A22 2 

w(k) w (k) 
en 
A 
Nes (k) = Air A22 ï (3.51) 
w(k) 

zodat Rr) = k (3.51) 


Karakteristieke versnellingstijd 


Voor (3.34) kunnen we nu met (3.36) en II(3.36) schrijven in het niet- 
relativistische geval (p = m v) 





BE ò dk kk òf 
BE a Br CW -eD EE den 
Ln an? de) kj of 
wt je ww 8 (k,, = he] (3.52) 
waarin ke, z (k.v)/v. 
Schrijven we 
dk 2 2 u 
[on WC) Bk), = WIN = (3.53) 


dan vinden we tenslotte voor de karakteristieke versnellingstijd voor niet= 


relativistische deeltjes uit (3.53) 


ed 


een 
A =l N. W £ 

gedafdd (En CDL Et & (3.54) 

k de aki Te Vte k, %pe we Vte k) 


In (3.53) zien we dus dat slechts die deeltjes die resoneren met aanwezige 
golven (ete) — k.v = 0) versneld worden (Cerenkov- of Landau resonantie). Uit 
(3.54) volgt dan dat versnelling van (niet-relativistische) deeltjes door 
Langmuir turbulentie zeer efficiënt ís. Aangezien belangrijke versnellings- 
mechanismen als bijv. Fermíversnelling (zie boven) en schokversnelling (zie 
beneden) slechts werken als de deeltjes aanvankelijk reeds een bepaalde 
“injektie-energie”" te boven gaan ‚valt te verwachten dat stochastische 
versnelling door plasmaturbulentie als “voorversneller" een belangrijke rol 


kan spelen. 


REGULIERE DEELTJESVERSNELLING 


b. Schokversnelling 


Bell (1978) en Axford (1976) hebben een versnellingsmechanisme in schokken 
voorgesteld. Het idee ís dat voldoende hoog-energetische deeltjes door de 
schok"diskontinuïteit” heen en weer kunnen reizen indien de botsingsfrequentie 
laag genoeg ís (lage dichtheid, hoge temperatuur). Indien er nu voldoende 
(elektro)magnetische fluktuaties aanwezig zijn aan beide zijden van de schok, 
zal de verdeling van hoog-energetische deeltjes mín of meer isotroop zijn in 
het referentiesysteem van het achtergrondgas. Maar het ongeschokte en 
geschokte gas bewegen zich aan weerszijden van de schok naar elkaar toe. Het 
komt er dus op neer dat enkele hoog-energetische deeltjes meerdere malen 
worden gereflekteerd als het ware tussen twee naar elkaar toe bewegende 
spiegels. Indien nu de ontsnappingskans van de versnelde deeltjes 
onafhankelijk van de energie ís ontstaat weer een machtsspektrum. Het ver= 
snellend mechanisme is . een eerste orde Fermimechanisme 
(dE/dt « LA ÎepeVve vt dat van de "tweede orde” ís) omdat we te maken hebben 
met de regelmatige beweging van twee naar elkaar toe bewegende spiegels waar- 
tussen alle spiegelende deeltjes energie winnen. Omdat de botsingsfrequentie 
laag is hebben we te maken met een zgn. botsingsloze schok waarbij het 
dissipatiemechanisme dat de schokovergang bewerkstelligt niet uit botsingen 
maar uit de een of andere vorm van plasmaturbulentie bestaat die de gerichte 


kinetische energie van de gasstroming vóór de schok omzet ín “verhitting” na 


I= 


de schok. Hoe deze plasmaturbulentie (dus de schokstruktuur) er nu precies uit 
ziet blijkt (gelukkig) niet zo belangrijk te zijn. 

We laten nu eerst ín ál. zien hoe magnetische fluktuaties (turbulentie) 
de verdeling hoogenergetische deeltjes isotroop zal houden, vervolgens ín 
b.2. hoe de gemiddelde snelheid van de hoog-energetische deeltjes er uitziet 
en tenslotte in 4.3. hoe een machtsspektrum tot stand komt. 


4.l. Spoedhoekdiffusie 


We beschouwen de interaktie tussen geladen deeljtes en elektromagnetische 


golven. De impulsverandering onder invloed van de golven is nu gegeven door: 


v x B 


d — _ é 
ge e= aôEta (4.1) 
We kunnen ÖE en ÔB weer ontwikkelen ín een Fourierreeks (= superpositie van 


vlakke golven) en kunnen (4.1) dan schrijven als: 


ä d° k dw ik.x(t) — tut 
ae =S Gpr kde (4.2) 


waarin 
Eko) Za Bk) to > 


B(k,‚w)} , (4.3) 


en x(t) weer de ongestoorde baan van het deeltje is. We kunnen F(k‚w) geheel 


in termen van E(k‚w) uitdrukken door gebruik te maken van de wet van Faraday: 
k x E(k,u) = = B(k, w) (4.4) 


Ek,u) za (Ek) tr * (k *ELW)} = 


Kev (v « E(k‚w)) 
nr BE en (4.5) 


(Bedenk dat A x (B xC) = (A . C) B = (A . B) C « ). Hieruit volgt dat qua 
grootte-orde de elektrische en magnetische kracht in een golf zich verhouden 
als: 

uE 


[M 
TEE xvj Ckv Mr 





de 


en dat de elektrische komponent verwaarloosd kan worden als w << kv of 
w= 0 (de zgn. magnetostatische benadering). 

De uiteindelijke vergelijking voor de evolutie van de deeltjesverdeling 
ziet er weer precies zo uit als (3,34), maar de uitdrukking voor D is nogal 
ingewikkeld, en we zullen hier alleen een speciaal geval beschouwen waar de 
elektrische kracht te verwaarlozen is. 

We beschouwen een storing in het magnetisch veld 


d° k dw ik.x =— Lut 
ÖB(x,t) = Í Tre Bk 0) e (4.7) 


en verwaarlozen het elektrische veld behorend bíj die storing. De bewegings 
vergelijking voor een deeltje luidt: 


dp _(v *Bo +v* 8) 


ee nn Rn 


waaruit direkt volgt dat de kinetische energie van het deeltje behouden ís 
tijdens de interaktie met de golf: 


2 dp 
dd pete. ger aep« (U *(B0 +65) = 0 (4.9) 


De ongestoorde baan ís nu een schroef (helix) rond het magnetisch veld. 


De bewegingsvergelijking is 


dp aq(v Xx Bo) 


de © c 
Er (4.10) 
dt ° tm 


Om deze op te lossen kiezen ge het magneetveld langs de z-as. Bog = Bg z. Uit 
(4.10) volgt direkt dat p . dE 
arbeid. De oplossing van het stelsel vergelijkingen (4.10) met als randvoor- 


= kr (Gp?) = 0. De Lorentzkracht verricht geen 


waarden x(0) = = (x05Y0»20) en v(O) = (t, ads Og, ve Ops vj) is dan 


Za 
gegeven door 


« 4 Bo 





a 


Ee, ko 
- + - + 
cos (Og Qt) x Te (Og - Qt) y v/2 








adh 


x= (zo - SF {ein 00 - sin (Oo - Q't)}) x + (4.11) 
a 
Vv. . de pn DN 
+ (yo + {cos Og = cos (Og = 2 0)}) y + (zo + vj) z 
a 


Voor de ongestoorde baan is de spoedhoek B, gedefinieerd door cos B = 
Pp « Bo/p Bo = p‚j/P» een behouden grootheid. . 

We beschouwen nu een storing B = ÖB xXx + o5, e die ís opgebouwd uit 
golven die zich langs het magneetveld ad nel dewez. k = k2. Uit (4.8) 


volgt met (4.9) dat de spoedhoek B onder invloed van ÖB verandert volgens: 


ôB êB 


d AB _ %* koe 08 _X 
er dek: { cos O Bj” Sin O5 } (4.12) 


p 

waar de hoek 0 = arctan (—) langs de ongestoorde baan moet worden berekend 
Ps 

(zie 4.11): 


%* 
O= O-Qt (4.13) 


We kunnen (4.12) schrijven in de vorm: 








ÖB êB 
dA * + -iO(t = +iO(t 
gend ee pe (4.14) 
waarin 
6B,(x(t),t) = 6B, (z(e),t) t êB,(x(t),e) (4.15) 
De berekening van ne gaat geheel analoog aan die van Stop in het geval 
van elektrostatische golven: 
tT + 1 gd 
<ABAB> = fadt' | ds 5 Bz {<8B, (£') ÖB,(t'+te)> x 
t 4 * 
=1(200-2 (2t'+s)) -iQ s 
xe + <öB, (t') ÖB_(t'+s)> e + 
* 
+iQ s 
+ <ôB (t') ÖB, (E'+5)D e + 
* 
+i(209 - 2 (2t'+s)) 
+ <8B_(t') 8B_(t'+s)> e } (4.16) 


De correlaties <ôB, CE!) ÖB,(E'+s)D etc. hangen alleen af van s (zie (3.8), wat 


betekent dat de eerste en de laatste term in het rechterlid van (4.16) weg- 
vallen bíj integratie over t' wanneer Tt > a, vanwege de term ai t' in 
de integrand. Invullen van de Fourier-ontwikkeling (4.7) levert dan na enig 
rekenwerk, en gebruik van 


Bike) Bj(k',w')> = „lim <B,(k,u) B; (k,w)> (27)* S(ktk') 5lwtw!) (W)TE, 


B ‚®% w) = B, (-k, —w) n, den „% w) B Tak, w)> = O0 vanwege cylinder-symmetrie 
Reek Bo: 
<|B, (‚| 2 
SABAB _ 1 o*2 dk dw 
ee nak je) (2m) 2 en 27 ölu-kv, +” ) + 
<|B_(k,w)| 2D 3 
NE ne id . 
+ 5 2n ôlurkv, „-Q >} (4.17) 
waarin: 
B,‚(k,w) = B (k‚w) £ 1 B Ck, U) ‚ (4.18) 


Vijg Bo” cos B. Merk op dat het verschijnen van de cyclotron=frequentie 
* 

Q = ymc 

afkomstig is uit de fase b= wt = k.x(t) = (w= kv) t van de golf zoals 





in de resonantie voorwaarde w- kv, „28 = 0 in dít geval niet 


gezien door het deeltje, maar uit de snelheidsafhankelijkheid van de Lorentz=- 
kracht, zoals blijkt uít (4.8) en (4.12). 


De Fokker-Planck vergelijking die: de evolutie van de distributiefunktie 
beschrijft ziet er ín dit geval als volgt uit: 





1 ô 
E- SIKB op (sta PD, 55) , (4,19) 
waarin: 
_  <ABAB> 
p= SER (4.20) 


Het verschijnen van factoren sin B in (4.19) ís een gevolg van het feit dat we 
niet met cartesische coördinaten ín de impulsruimte werken, maar met p;, 
B en 0. 
Vergelijking (4.19) beschrijft een diffusie van spoedhoek van de 
deeltjes, ook wel “"spoedhoek verstrooiing” (pitch angle scattering) genoemd. 
Het effekt ervan is dat een verdeling f(p,B) die anisotroop is, d.w.z. waarin 


er een voorkeursrichting voor de komponent van de impuls p langs het magneet 


Et De ne Ee ek ld ie did 
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veld is, onder ne van die diffusie zal streven naar een isotrope 







verdeling GERELN 5e En = 0, (zie figuur). 










n veaukromme 5 


SOp 0) zest 





anusotfroofp Ísofroop 
Dit ís het eenvoudigst als volgt in te zien: 


definië®r: 
gn ke 
fp) = k | dB sin B £(p,B) 
o 
7 em 
I= k J dB sin B {£(p,B) - E(p)}2. 
Le) 


Dan volgt uit (4.19) 


ÒI ÒÉ 
halle j dB sin B Dg Ga) < 0. 


Voor iedere anisotrope verdeling hebben we I > 0, en de spoedhoek verstrooiing 


reduceert die anisotropie. 


4.l.1l. Dit proces ís bijv. van belang bij de voortplanting van snelle deeltjes 


door de zonnewind. Men neemt daar fluktuaties ín het magneetveld waar, en deze 


fluktuaties verstrooien de snelle deeltjes, zodat hun verdeling (bijna) 
isotroop wordt. Dit beperkt de voortplantingssnelheid van die deeltjes langs 


het magneetveld: 


pe 

<v,,> = k | sin B dB v cos B f(p,B) = 
HH at 
+1 


=k Ĳ, d(cos B) v cos Bf (p,B) = 
+1 


= bk f d (cos B) „Leef tE LmD, (4.21) 
1 cos B 


Hieruit blijkt direkt dat voor een isotrope verdeling Cv, jp = 0, en dat ín 
het algemeen de gemiddelde voortplantingssnelheid evenredig is met 


=27- 


ÒË 


òf 
Boog p "sin Pie 


4.1.2. Voor links en rechts circulaire gepolariseerde Alfvêén golven 


Gu = + 1 respektievelijk) geldt: 
y 
B, (k, w) = k B, (k) (z = 1 2n ôlw=- kv) 
B, (k,w) = % Br (k) ûz + 15) 2n ô(w= kv‚) (4.22) 


* 
en kunnen we uit (4.17) afleiden voor het geval kv, << Q +] kv | 


„ak Bo? „ « klB,(o) |? 
Dg =, 8 ET ie d +78 ET 8 (4,23) 
vy B. 


4.2. Diffusie van snelle deelt jes langs het magneetveld 


Uit (3.17) volgt, dat een niet relativistisch deeltje, dat door een 

geioniseerd gas bestaande uít elektronen en protonen loopt met een snelheid v 

die veel groter ís dan de thermische snelheid van de elektronen en protonen, 

met de deelt jes in het gas botst met botsingsfrequentie 
ee 2 Bit a 

(E = 5 mv + mc“; Ve De IlnA (ár N) (zie II.1.1) 


ws ì 
ik. 
vw) Va e= Je (4.24) 
kT 5 
waar v, = Er de thermische snelheid is. Dit betekent dat deeltjes met 


k_ T k_ T 
snelheid wv >Ò Gerek: Ey vrijwel ongehinderd door botsingen door het 
plasma lopen. In een plasma met een magneetveld blijft deze conclusie 
ongewijzigd voor beweging langs het veld. 

We hebben echter al gezien dat kleine fluktuaties in het magneetveld ín 
staat zijn de spoedhoek 6 van een deeltje te veranderen (zie (4.19). Als 
zodanig spelen ze dezelfde rol als botsingen, en men zou kunnen verwachten dat 
daardoor de deeltjesvoortplanting langs een magneetveld Bo = Bg ë onder 
invloed van spoedhoekverstrooiing ook een diffusieproces ís met diffusie- 


coëfficiënt: 


2 


B 


wol 


1 
AT ke Baakdn (4.25) 


=28- 


_ SABA 
T 


waar v de rol van de botsingsfrequentie speelt. We zullen deze 


=D 
B 
relatie formeel afleiden. We voeren ín u = cos B, en nemen aan dat de 
verdeling van energetische deeltjes een gradiënt 2 langs het magneetveld 


vertoont. In dat geval wordt vergelijking (4.19) 


(4.26) 


Bett VPB” Bu tZ You 


<ABABD 
T 


ek DÉ _ À de dE) 


waarin v = ‚ We nemen aan dat het verstrooíiiíingsproces dominant is, 


v > », en ontwikkelen f als: 





mf Eje 0(e?) (4.27) 
fl 1 =1 
waar 7 O(e) en € aak {<1 ‚ In order e volgt dan: 
o. 
d A - fe) f òf, 
BE dp v on 2) = 0 > Dr hd 0 » (4.28) 


wat betekent dat de verdelingsfunktie isotroop is tot orde ed. Gebruiken we 
dit resultaat in (4.26) dan blijkt na middeling over u: 


+1 DE, BE 
5 | sr W-g=d. (4.29) 


In orde eÔ tevert (4.26) 


nt 2 fi 
ba “D Cn nk. on 5) . | (4.30) 


Na integratie over p van =l tot u volgt: 


2 ÒÉ 2 Of 1 
pn 1 a 
bv Le EE war (4.31) 





met als oplossing: 


f tv 
Fi Eil) ee | dr er a  « (4.32) 
zi v(u') ôz 
In orde e levert (4.26): 
Òf Of 2 fz 
- RS en 
rrd vi 5 ar au 7 rn Je (4.33) 
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Invullen van (4.32) en middelen over pu (aangegeven door =) levert: 


| 


Ò ‘1 v te 
amen and meme ' an ee EE Hi 
seb 8 du vu A du rj 0 (4.34) 


Door partiële integratie kunnen we (4.34) schrijven als: 


EN of 


ze 0 Se 
met 
7 dei 
D,/ zk Í du v 2D) (4.36) 


Uit (4.29) en (4.35) volgt dan dat de “transportvergelijking" voor 
F=f + Éj 
o 


òf _ ô rd 
ôt “ dz P// oz AE 


waar we f, hebben vervangen door E,‚ wat een fout van orde e? introduceert, die 


niet van belang ís omdat Ef gedefinieerd is tot orde €. 


Wanneer v(u) = v, onafhankelijk van u is, kunnen we D,, precies 


B 


uitrekenen en volgt 


2 
(4.38) 


<|s 


wl 
» 


Dip 


in overeenstemming met (4.25) uit de meer heuristische beschouwing aan het 
begin. 

De konklusie is dat energetische deeltjes, die op grond van botsingen met 
thermische elektronen en protonen niet kunnen worden gestoord in de beweging 
langs het magneetveld, toch niet vrij kunnen bewegen als er fluktuaties ín het 
magneetveld zijn. De effektieve botsingsfrequentie ten gevolge van fluktuaties 


met amplitude êB volgt uit (4.17) en is qua grootte-orde: 


7E d | (439) 


4,3, Versnelling door stationaire quasi=evenwijdige schokken 
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4.3.1. Parallelle schokken 


We bekijken het geval van schokken die zich evenwijdig aan het achtergronds- 
magneetveld voortplanten. Het in de inleiding geschetste schokversnellings- 
mechanisme werkt onder twee voorwaarden: 

l. Energetische deeltjes zien de schok als een “passieve” diskontinuïteit én 
2. Vóór en na de schok zijn er voldoende magnetische fluktuaties die de snelle 
deeltjes praktisch isotroop maken in een referentiesysteem dat in snelheid 
niet wezenlijk verschilt van het betreffende rustsysteem van het geschokte dan 
wel het ongeschokte gas. Het is mogelijk aan te tonen dat de snelle deeltjes 
die stroomopwaarts ontsnappen door hun anisotrope verdeling zelf fluktuaties 
opwekken ín de vorm van Alfvéngolven die van de schok af gericht zijn (De 
groeisnelheid van een dergelijke instabiliteit kan men berekenen met de 
uitdrukking voor de- groeisnelheid in I1(3.40).) De verdeling van snelle 
deeltjes stroomopwaarts wordt dan praktisch isotroop in het golfraam dat zich 
met snelheid w‚-v, ten opzichte van de schok beweegt (zie figuur beneden). 
Echter ook stroomopwaarts dienen voldoende magnetische fluktuaties te bestaan 
zodat snelle deeltjes een grote kans hebben gereflekteerd te worden. We 
onderstellen nu voor het gemak dat zowel vóór als achter de schok voldoende 
turbulentie aanwezig is om de snelle deeltjesverdelingen isotroop te houden 
zodat we 4.1 en 4.2 kunnen gebruiken. Bekijk de figuur waarin de stroming is 
aangegeven voor een parallelle schok (d.w.z. met schoknormaal parallel aan het 
magnetisch veld) ín het referentiesysteem van de schok (zodat de schok 


stilstaat). 


W nn Za wa @) 





Schele 


Grootheden in het ongeschokte gas geven we aan met index 1, die in het 
geschokte gas met index 2. Het ongeschokte gas beweegt zich naar de schok toe 


met relatieve snelheid wis het geschokte gas beweegt zich van de schok af met 
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snelheid w, < Wi (Merk op dat het in het algemeen mogelijk fs uit het behoud 
van massaflux, impulsflux en energieflux door de (stationair veronderstelde) 
schok heen relaties tussen de grootheden vóór en nà de schok op te stellen 
indien de materiaaleigenschappen van het plasma bekend zijn; biíjv. Rankine- 
Hugoniot relaties voor een ideaal gas met adiabatische index y). Het plasma 
voert magnetische fluktuaties met zich mee, die deeltjes verstrooien. We nemen 
aan dat deze fluktuaties stil staan ten opzichte van het plasma. In het 
stelsel waarin het plasma in rust is, hebben deze fluktuaties geen elektrisch 
veld (6E 2= ze Xx 8B) ), en ís de botsing elastisch, d.w.z. alleen de spoed- 
hoek 8 tussen het magnetisch veld en de snelheidsvector van het deeltje ver- 
andert, maar niet de energie. We noemen de plasma ruststelsels vóór en achter 
de schok S; respectievelijk Sj. We beschouwen nu een deeltje dat van achter de 
schok onder een hoek B, de schok kruist naar gebied 1. We nemen aan dat het 
deeltje een impuls p‚ en energie E ín het stelsel S) heeft. In het stelsel S, 


heeft het dan een energie (zie Lorentztransformaties II.3.8). 
E, = YE _ Be c) = ‚A (E-(u,-v)p cos B, (4.40a) 


en een impuls 


Pi “ Y, TB Be) = velp — vm(wj-w)z (4.40b) 
waarin An w‚-w)/e, le (1-8, en y= (L-v2/c2) | met v de deelt jes- 
snelheid. 


Als het deeltje nu in Sj elastisch wordt teruggekaatst zonder dat zijn 
energie dus in S, verandert, en de schok kruist onder een hoek B dan is 


zijn energie in het stelsel 55: 


En 1. F4 + (w, = w‚) py cos B) (4.41) 


We nemen aan dat | =v2 | << ec en dat p >> |ym w‚-w)l zodat tot in orde 
(ww) /v << 1 de energieverandering is 


AE = E, „EE = (ww) p(cos B =cos BĲ) (4.42) 


We nemen aan dat de spoedhoekverstrooiing door de magnetische fluktuaties 


sterk genoeg is om de deeltjesverdeling vóór en achter de schok bijna isotroop 


Ming PE de OT Tt EV WE RAE 7 Wi SE Dd 


En 


te houden in het stelsel Sj en respectievelijk Sj. Het aantal deeltjes dat van 
gebied 2 naar gebied 1 kan oversteken onder een hoek tussen B, en 8, + dpj per 
tijdseenheid en per eenheidsoppervlak van de schok is dan evenredig met 27 


cos Bj sin Bj dBi, en de gemiddelde waarde van cos fj over alle hoeken is dus: 





T « 
2x cos? Bj sin Bj dB) 
cos B = =-à (4.43) 
Á 2ncos B, sin B; dB, 
Er 
Evenzo volgt cos Bj = + 2. De over de hoeken gemiddelde energietoename is 
dus: 
dend (wi-wa) v ® 
a B Me B Aant A (4.44) 


De deeltjes winnen dus energie geheel analoog aan het geval van een deeltje 
dat opgesloten zit tussen twee met snelheid wj = wz naar elkaar toe bewegende 
spiegels. De enige funktie van de schok is, dat hij de gewenste diskontinuï- 
teit in de plasmastroming veroorzaakt. Hieruit volgt al meteen, dat; alléén 
deeltjes met een vrije weglengte voor spoedhoekverstrooiing \ = GÉ) >> D, 

waarbij D de dikte van de schok is, versneld worden. Is het omgekeerde het 
geval, dan “zien” de deeltjes geen discontinuïteit ín de snelheid van het 
plasma (en de veldfluktuaties). Wanneer per tijdseenheid de schok v maal wordt 


overgestoken, dan neemt de energie van een deeltje toe als: 


== (wi = wa) v 
dE > 4 
re Miamnne” aamendn (4.45) 


Natuurlijk blijven de versnelde deeltjes niet altijd bij de schok, maar ont= 
snappen na verloop van tijd en worden door het plasma meegesleept naar z 
= to, Vergelijking met (1.11) laat zien dat (4.45) op te vatten ís als 


reguliere (i.e. eerste orde in (wj = wo)) Fermiversnelling. 


Het spektrum: De voortplanting van de snelle deeltjes voor en achter de schok 
vinden we nu door het effekt van “meesleping” (konvektie) door het 


achtergrondgas aan (4.37) toe te voegen: 





B) 4 8 (NE) =D), NE) = 0 (4.46) 
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Hierin is NC) dE het aantal deeltjes met energie tussen E en E + dE en 


wol 


D,, ad a de diffusiecoëfficiënt langs het magneetveld zoals berekend in 
(4.38). Dé fysische verklaring van (4.46) is de volgende: spoedhoekver- 
strooiing door de magnetische fluktuaties tracht de verdeling van de snelle 
deeltjes isotroop te houden ín het stelsel dat met de fluktuaties (en het 
plasma) meebeweegt. Als dat volledig zou lukken zouden de deeltjes in dat 
stelsel geen gemiddelde snelheid langs het veld hebben, en de deeltjesflux ten 
opzichte van een stilstaande (t.o.v. de schok) waarnemer was gelijk aan 
Fo = w NE) (j = 1,2). Echter, een gradiënt in de deeltjesdichtheid intro-. 


duceert een kleine anisotropie ín het meebewegend stelsel (zie (4.31) met een 








daarbij horende flux ten gevolge van diffusie F, == D,/ En ‚. De totale 
deeltjesflux is dus gelijk aan F(E) = Fo + F, = WN(E) -— D,, eg ) ‚ en (4.46) 


beschrijft deeltjesbehoud. Als er geen tijdsafhankelijkheid ís volgt uit 
(4.46) 


wN(E) — D LE) = constant , (4.47) 


met als algemene oplossing: 
z 
N(E,z) = A +B exp [{ dz w/D,,1 (4.48) 
o 
Voor de schok (z < 0) wordt dan de oplossing: 
. z 8 
N(E,z) = N(E,-e) + {N(E,O) = N(E,-e)} exp [{ dz wal, ,l (4.49) 
o 


en de flux ís gegeven door: 


Achter de schok is de enig fysísch acceptabele oplossing (z *» +e): 


N(E,z) = N(E,0) 
(4.51) 
F3(E) = wz N(E,o) 
Vóór de schok ís de gemiddelde deeltjesflux naar de schok toe gericht, en alle 
deeltjes die in dat gebied terecht komen zullen de schok weer een keer 
bereiken. Deeltjes kunnen dus alleen aan het versnellingsproces ontsnappen 


achter de schok, waar de gemiddelde flux van de schok af gericht is. In goede 


me, nn 


Wz 
benadering S= << 1) is de flux van deeltjes over de schok vanuit gebied 2 


‘naar 1 (bedenk dat N(E) dE = ár £ p? dp) 


& 
r.=d f 2m v cos BN(E) sin Bdp=-LvNE) . (4.52) 
+ án on 4 
De ontsnappingskans n is nu gegeven door: 
ee en << 1 (4.53) 
n bed mmm ne __ . . 
F + v 


. n 
De kans „ve n maal de schok over te steken is dus P, = (l-n) =1=-an= 





1 = y__» aangenomen dat elke oversteek onafhankelijk van de vorige ís, 
zodat (nn << 1) 


In P, = an | (4.54) 


Na n oversteken ís de energie van een extreem relativistische deeltje (v 


= C) met oorspronkelijke energie E, volgens (4.97) toegenomen tot 
n 
E. = E, (1+a) (4.55) 
zodat (a << 1) 
An (E/E) = nd. : (4.56) 


Aangezien ín een stationaire situatie het aantal deeltjes NE) met energie 


E > E‚ evenredig is met P_, dus 
An NIE ) = Kk + An P 
n n 


vinden we tenslotte 


En E n/a 
Hiele (4-57) 
(4.56)  ° 


zodat het differentiële energiespektrum de vorm heeft 


2wg + wi 
—(1l+n/a) 7 Wi = W2 
NE) = UE , E-) = (E) (4.58) 


0 o 
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Zowel voor supernovaresten als voor extragalaktische radtiobronnen (hot spots) 

lijkt dit schokversnellingsproces hese Voor sterke schokken kan men 
— TEL 

afleiden dat — "ty a ’ 4 , waar Y 3 de adiabatische index is. In dat 


geval is de Me 2 (vergelijk de waargenomen waarde 2.5 = 3). 


4,3.1l.1. Injektieënergie 


Ook indien Coulombverliezen niet belangrijk zijn (zie 1.3.) is er een 
“injektieënergie” nodig indien de magnetische fluktuaties uit Alfvéngolven 
bestaan. In dat geval treedt verstrooiing van deeltjes namelijk slechts op als 


een van de resonantiekondittes (zie (4.,17)) 


* 
w=- k,/vyy Q = 0 (4.59) 
is voldaan (mêt de juiste polarisatie!). Voor Alfvéngolven (w= kj, vi) volgt 
dan Cv, << v‚y ondersteld) 
* 
v‚j| = v|cos B|= v [& (4.60) 
[v‚/ A Vo , 


zodat voor deeltjes van soort j 








En 
Zes hell ae? 25 m0? | (4.61) 
j Va w j 1 


es =l < 
(met Bd = q3 B, me) ) omdat voor Alfvêngolven w w,/2: Voor elektronen 


en protonen levert dit een minimale injektieimpuls. 


Pe “ P4 “ E, v/e? & 2m, Ve (4.62) 


Vooral voor elektronen lijkt dus voorversnelling gewenst. 


De Konstante elektrische velden 


Versnelling van deeltjes geschiedt in een plasma uiteindelijk door elektrísche 
velden (de Lorentzkracht verandert de deeltjesenergie niet). Konstante 


elektrische velden zijn zeer efficiënt ín het versnellen van deeltjes mits het 
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elektrisch veld gealigneerd ís met het magnetische veld. Voor elektrische 
velden loodrecht op het magnetische veld treedt nauwelijks versnelling op 
zolang clE | < vil, Omdat het elektrische veld veelal door bewegingen ontstaat is 
echter |E | ha 8 _X B | <IEl. Deze overwegingen beperken het voorkomen van 
versnelling door konstante elektrische velden ín astrofysische situaties 
aanzienlijk. Een voor versnelling belangrijke situatie doet zich voor indien 
rekonnektie van magnetisch veld optreedt. Plaatselijk waar het magnetische 
veld verdwijnt ("magnetisch neutraal gebied") of van richting verandert 
(“stroomlaag”) kunnen konstante elekrische velden (als gevolg van 
fluxannihilatie (v xE= zer ÒB/òt) of als gevolg van 
potentiaalverschillen díe elders door bewegingen ontstaan (E = B XB) en in 
het gebied van rekonnektie worden “kortgesloten") optreden. We bekijken nu de 


situatie E//B. 


5.l. Runaway versnelling 


In een thermisch plasma van dichtheid n en temperatuur T zal een konstant 
elektrisch veld (dus //B genomen) een elektrische stroom oproepen. Gemiddeld 
zal de energie die door een “thermisch” elektron na een Coulombbotsing uit het 


elektrisch veld wordt opgenomen bij de volgende botsing weer worden afgestaan: 


dv eE 
o= TE = <= el VP ; (5.1) 


Er ontstaat dus een gemiddelde driftsnelheid van de elektronen (en een veel 


kleinere van de zwaaardere ionen) die uit (5.1) volgt: 


® 
ls | 


Vp Sm. (5.2) 
c 
w An A 
waarin de botsingsfrequentie Vn er Ci Kie (vgl. II 1.1). (5.3) 


Aangezien de deeltjesverdeling een Maxwellse vorm heeft zijn er echter ook 
enkele snellopers met v > Var Voor deze deeltjes geldt (zie (1.17)) 

dv eE 
e= ma _ vv) yv (5.4) 


met 
v.v, +). (5.5) 


== 


Uit (5.4) volgt nu dat deeltjes met 


ve. Vie nm 5 e w B An A 5% 
vi) = (5.6) 


e E E 
voortdurend zullen worden versneld, of althans niet op grond van Coulomb- 
botsingen voldoende wrijving voor een stariótaire deifesnetheid ondervinden. 
Zodra het elektrisch veld zo groot wordt dat de wegloopsnelheid (5.6) zakt tot 
de thermische elektronensnelheid neemt de versnelling “ernstige” vormen en zal 
het gehele plasma worden versneld. Dit gebeurt dus bij de kritische veld- 


sterkte (ongeveer tweemaal de zgn. Dreicerveldsterkte) 
ee 2 
E. e ko mÂ (5.7) 


Met (5.7) kunnen we nu de kritische snelheid (5.6) waarboven “runaway” 
optreedt herschrijven als 


E& 
vS) v (5.8) 


c te 
Omdat de deelt jessnelheid altijd kleiner dan de lichtsnelheid is, voglt uit 
(5.8) dat runaway slechts kan optreden als 
ko er 


E 
Eme? (5.9) 


5.2. Dubbellagen 


Een belangrijk fenomeen vormen de nog slecht begrepen elektrostatische dubbel- 
lagen; deze zijn voor te stellen als koherente elektrostatische golven van 
grote amplitude met sterke elektrische velden ten gevolge van ladings- 
scheiding. Mogelijk spelen zij een belangrijke rol bij de versnelling van 


deeltjes in de astrofysica. 
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